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Abstract The observation of the interstellar 1.809 MeV decay-line of radioactive 
26 Al by the imaging gamma-ray telescope COMPTEL have let to the 
conclusion, that massive stars and their subsequent core-collapse su- 
pernovae are the dominant sources of the interstellar 26 Al abundance. 
Massive stars are known to affect the surrounding interstellar medium 
by their energetic stellar winds and by the emission of ionising radia- 
tion. We present a population synthesis model allowing the correlated 
investigation of the gamma-ray emission characteristics with integrated 
matter, kinetic energy and extreme ultra-violet radiation emission of as- 
sociations of massive stars. We study the time evolution of the various 
observables. In addition, we discuss systematic as well as statistical un- 
certainties affecting the model. Beside uncertainties in the input stellar 
physics such as stellar rotation, mass loss rates or internal mixing mod- 
ifications due to a unknown binary component may lead to significant 
uncertainties. 



1. Introduction 

During recent years spatially resolved observations of decaying ra- 
dioactive isotopes spread in the interstellar medium of the Milky Way 
by their characteristic 7-ray emission became feasible. The Compton 
telescope COMPTEL [1] aboard NASA's Compton Gamma-Ray Obser- 
vatory surveyed the sky in the MeV regime during it's 9 year mission 
time. One of the key results of this mission was the generation of the 
first all sky image in the 1.809 MeV 7-ray line from radioactive 26 Al 
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[2, 3]. The COMPTEL observations confirm a integral 26 Al content of 
the Milky Way of 2 to 3M Q . In principle 26 Al can be produced by 
various nucleosynthesis sites. Since it's first detection by HEAO-C [4] 
non-explosive sites such as AGB stars [5] and Wolf-Rayet stars [6] as 
well as explosive environments (novae [7] and core-collapse supernovae 
[8, 9]) are discussed as possible source candidates for interstellar 26 Al. 
The reconstructed 7-ray intensity distribution correlates best with trac- 
ers of massive stars such as thermal free-free emission from the ionised 
medium [11]. Due to its long lifetime freshly released 26 Al travels quite 
some distances until its final decay. Therefore the 1.809 MeV line emis- 
sion is rather diffusive and it was not possible to observe isolated can- 
didate sources. 7 2 Velorum, which is the nearest Wolf-Rayet star, is 
the only possible source for which significant upper limits on the 1.8 
MeV flux could be extracted so far [12]. Because of the close correla- 
tion between the galactic distribution of massive stars and the observed 
1.8 MeV intensity pattern, aggregations of massive stars such as young 
open clusters and OB associations seem to be suitable laboratories to 
study these candidate sources. In particular, the by combined analysis 
of the measured 1.8 MeV intensity together with observables in other 
wavelength bands one may gain insight in the physical processes in- 
volved. Massive stars not only eject freshly synthesised material into 
the surrounding ISM but also impart a huge amount of kinetic energy 
by stellar winds and their subsequent core-collapse supernovae. Fur- 
thermore they affect the state of the interstellar medium by emitting 
a large fraction of their electromagnetic radiation in the photoionizing 
extreme ultra-violet regime. Therefore correlations between the 7-ray 
line emission of 26 Al and observables being related to one or the other 
phenomenon of massive stars are expected. 

We make use of a detailed population synthesis model in combination 
with an 1-dimensional thin shell expansion model of evolving superbub- 
bles to study the time evolution of the different emission parameters 
and their resultant observables. In the following sections we present the 
model and discuss the predictions and their uncertainties in detail. Fi- 
nally, we apply the model to OB associations in the Cygnus region and 
discuss the results for the expected 26 Al content of this area. 

2. Modelling the OB star Emission 

2.1. Release of radioactive Isotopes 

26 Al is synthesised via proton capture onto 25 Mg and therefore is a 
secondary product of the CNO-cycle in massive stars. The efficient pro- 
duction of 26 Al requires temperatures of at least 3.5 x 10 7 K, which are 
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reached by stars initially more massive than 25M Q . The production rate 
increases as the core temperature approaches values of ~ 10 s K and de- 
pends strongly on the initial metallicity of the stellar material [6]. During 
later burning stages 26 Al is produced by proton captures of secondary 
protons onto freshly synthesised 25 Mg. The main destruction channels 
at temperatures below 220 keV are the /3 + -decay and subsequently pro- 
ceeding reactions with 26 Al as input whereas for the late burning phases 
the direct decay of the isomeric form of 26 Al which is 220 keV above the 
ground-state destroys the remaining 26A1 very efficiently The internal 
mixing processes bring freshly synthesised 26 Al into areas of the massive 
star which will be expelled by the strong stellar winds during the Of- and 
WR-phases. Meynet et al. [6] computed a stellar evolution grid with 
detailed nucleosynthesis of 26 Al. Their yields could be sufficiently well 
described by a power-law fit, which is used in our population synthesis 
model to describe the 26 Al release of Wolf-Rayet stars. The uncertain- 
ties in the description of the mixing processes, stellar rotation and mass 
loss give raise to uncertainties in the expected yields of factors up to 3. 
Beside Wolf-Rayet stars core-collapse supernovae are demonstrated to 
be efficient sources of interstellar 26 Al [8, 9, 10]. In these events 26 Al 
is produced by hydrostatic and explosive Ne-burning as well as the v- 
process. Whereas for the mass range below 3OM detailed explosive 
nucleosynthesis models exist, which cover a wide range of input physics 
including different mixing schemes and stellar rotation, for type Ib/c 
supernovae this aspect is rather unexplored. The typical type II yields 
are fit reasonably well by a power-law. In contrast, the convergence of 
the core-masses for the most massive stars lead to the assumption of a 
non-mass-dependent 26A1 yield, which is supported by the findings of 
Woosley et al. [9]. As in the case of the WR contribution the theo- 
retical uncertainties give raise to an inaccuracy of a factor 2 to 3. In 
addition, core-collapse supernovae release a large amount of radioactive 
60 Fe, which migth be observed by the 7-ray lines of its daughter nucleus 
60 Co at 1.173 and 1.332 MeV. Due to its solely production by supernovae 
and the longer lifetime, which reduces the intensity of the expected 7- 
rays, the Ironlines are still undetected so far. Nevertheless, the ejection 
of 60 Fe is included in our population synthesis model. 
Figure 1 shows the time evolution of interstellar content of 26 Al for an 
OB association assuming an instantaneous starburst with Salpeter-IMF. 
The resulting time profile is convolution of the mass spectrum and the 
mass- and time-dependent ejection rates of 26 Al taking the radioactive 
decay into account. The profile shows the typical two-peaked struc- 
ture, where the first peak is due to the WR contribution and the second 
originates from the SN activity. Usually OB associations have limited 
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population statistics, leaving the observer with a poor sampling of the 
initial mass spectrum. A Monte Carlo version of our population syn- 
thesis code allows the estimation of the uncertainties introduced by this 
effect. Our finding is, that for populations being initially richer than 100 
stars more massive than 8M Q , the resulting time profiles are dominated 
by theoretical uncertainties instead of statistical effects. 
An additional source of uncertainty is a possible contribution to the in- 
terstellar 26 Al content by peculiar massive close binaries. Langer et al. 
[13] have shown that some of these systems may give raise to a large 
overproduction of 26 Al in the supernova explosion of secondary. Their 
contribution may be enhanced by factors up to 1000. Up to now, it is 
not possible to really incorporate these effects in a population synthe- 
sis model due to a lack of a detailed investigation of the appropriate 
parameter space. Therefore we have to leave this as open question. 



A well known phenomenon of massive stars is their very strong stellar 
wind imparting a fairly large amount of kinetic energy to the ISM. In 
addition, their subsequent supernova explosions release on a very short 
time-scale even larger amounts of energy to the ISM. In OB associations 
these effects can be expected to overlap and serve as an energy source 
for blowing large-scale gas structures, known as superbubbles, into the 
surrounding ISM. Our population synthesis model predicts the flux of 
kinetic energy as well as matter due to an underlying massive star popu- 
lation. The stellar wind part is modelled using the semi-empirical mass 
loss rates from the Geneva stellar evolution models [14] in combination 
with a semi-analytical wind velocity formula from Prinja and colleagues 
[15, 16]. The energy released by a single supernova is assumed to be 
10 51 erg per event. By using a similar scheme as for the nucleosynthesis 
part our population synthesis model predicts the integral flux of energy 
and matter from an association. 

The integral wind power as well as the flux of matter can only be ob- 
served by their impact on the surrounding ISM, therefore we used the re- 
sults from our population synthesis model as input for an 1-dimensional, 
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numerical model of an expanding supershell. The model is based on the 
thin shell approximation and incorporates radiative cooling as well as a 
parameterised description of evaporation from the inner shell or cloudlets 
which passed the shock and made it into the hot bubble medium. Our 
description of a possible evaporative poisoning as it was labelled by Shull 
& Saken [17] aims on a study of the relevance of this process and uses two 
parameters, which are the transmission efficiency e por os and the evapo- 
ration time-scale T eva p. Depending on the ambient density and pressure 
evaporation becomes critical for quite low transmission efficiencies of the 
order of few thousands relative to the swept up mass sitting in the thin 
shell. This is especially true if the evaporation time-scale is short. If 
the bubble medium reaches a critical density at sufficiently low tem- 
peratures of some 10 6 K the interior energy of the bubble is converted 
into radiation and the bubble medium is cold down to some 10 4 K rather 
immediately. If the evaporation time-scale is long compared with the dy- 
namic time-scales then the transmission of some fractions of the ambient 
medium reduces the mass of the shell at a given time and the bubble 
expansion stays faster than in cases where no material is transmitted 
through the shell. In this cases the bubbles become larger. Figure 2 
shows the time evolution of the ratio of the observable kinetic energy 
of an expanding supershell relative to the kinetic energy from source. 
In this case intermediate values have been chosen for the evaporation 
parameters. 

Figure 2 Ratio of observ- 
able to injected kinetic en- 
ergy of an expanding su- 
perbubble around an OB 
association. The surround- 
ing density was chosen to 
100 cm-3 whereas the evap- 
oration parameters were 
o 10 20 30 40 set to T poros = 0.005 and 

t[Myrl r cvap = 0.5Myr. 

During the early phase the ratio approaches the canonical value of 20- 
25% rather quickly, whereas after approximately 15 Myr the ratio drops 
significantly and approaches values of 5% or less in the late phase. This 
behaviour depends strongly on the chosen parameter set for the envi- 
ronment as well as evaporation process. If evaporation is efficiently sup- 
pressed due to magnetic fields or some other processes the ratio stays 
near 20% over the whole evolution. But, the critical parameter range for 
efficient cooling and subsequently stalling the expansion are limited by 
very modest values, which might be typical for interstellar conditions. 
The combination of these gas dynamic effects with the emission charac- 
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teristics for the release of freshly produced radioactive isotopes such as 
26 Al and 60 Fe allows to predict the spatial intensity distribution of the 
expected 7-ray emission. This is of special importance in the physical 
understanding why the observed 1.809 MeV intensity distribution cor- 
relates so well with strongly blurred tracers. Beside the uncertainties 
due to the evaporation process and therefore the efficiency of radiative 
cooling in the early stages, the results of the population synthesis are 
strongly affected by the uncertainties in the stellar wind models. The 
Geneva stellar evolution which have been used for consistency reasons 
apply an enhancement factor of 2 to the mass loss rates during the main- 
sequence and WR-phase. Therefore the expected wind power is roughly 
a factor of 2 higher than for typical observed stars. In contrast, an 
underestimation of the wind velocities of only 30% will already restore 
the wind power in our model. We therefore estimate our results to be 
reliable within a factor of 2 or so. 

2.3. Ionising Radiation 

As mentioned in the introduction the observed 1.809 MeV intensity 
pattern correlates very nicely with the galactic free-free emission ob- 
served with COBE DMR at 53 GHz [11, 18]. The interpretation of this 
close correlation points to massive, hot stars being the dominant sources 
of interstellar 26 Al. The free-free emission originates from the ionised 
portions of the interstellar medium which is found in compact HII regions 
as well as in from of the diffuse ionised medium. It is strongly believed 
that in both cases the ionising extreme ultra-violet radiation from young, 
massive stars is the dominant source of keeping the medium ionised [19]. 
In the framework of our population synthesis model we therefore incor- 
porated the emission of Lyman continuum photons by means of a simple 
fit-function in dependence of the initial mass of the emitting star. 



This function is fitted to the Lyman continuum photon fluxes from de- 
tailed stellar atmosphere calculations [20]. For main-sequence stars the 
error is less then 5% by using this fit instead of the appropriate model 
values for the given initial mass. Nevertheless, the error of the resulting 
time-profile of the Lyman continuum emission of an OB association is 
considerably larger due to disregarding the effects of stellar evolution. 
However, a comparison with an alternative EUV emission model based 
on direct interpolation of the appropriate fluxes from stellar atmosphere 
models reveals an overestimation of the cumulative Lyc flux in our as- 
sociation model of less than 30% which is of the order as theoretical 
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uncertainties of the underlying stellar models. 



Figure 3 Evolution of the 
correlation of the Lyman 
continuum flux vs. the de- 
cay rate of interstellar 26 Al 
and 60 Fe, respectively 
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Figure 3 shows the time-evolution of the integral Lyc flux versus the 
decay rate of interstellar 26 Al and 60 Fe, respectively, for the same as- 
sociation model discussed earlier. The plotted trajectories start with 
low decay rates at high EUV luminosities and evolve to considerable de- 
cay rates at vanishing Lyc fluxes. In principle, this behaviour could be 
exploit to construct a very sensitive age indicator for the underlying pop- 
ulation. The ratio of the observed 1.809 MeV flux from an astronomical 
source population and the respective Lyc flux, which in principle could 
be extracted from the observed free-free intensity by using a proper ion- 
isation model, should be independent of the distance to the population. 
After an additional normalisation of this ratio relative to the Lyc flux 
of an 07V star one gets the 07V equivalent yield for an radioactive 
isotope for a given population [21]. Indeed, this quantity can be used as 
sensitive age indicator as shown in figure 4 for the case of 26 Al. 



Figure 4 Time evolution 
of the 07V equivalent yield 
for 26 Al in comparison to 
the age ranges for other age 
indicators in astronomical 
populations. 
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Especially in the regime between 3 and 10 Myr the 07V equivalent yield 
shows a strong increase over 3 orders of magnitude. This is due to the 
strong decrease of the number of O stars because of supernovae and the 
two emission peaks in the time profile of the interstellar 26 Al content (cf. 
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Fig. 1). 

Additionally, the released Lyc flux can be used as input to an ionisa- 
tion model to calculate a prediction map of expected free-free intensity 
of an area under investigation. This could in turn be compared to the 
observed one and therefore may give an additional constrain on the real 
population. 

3. Statistics 

As already discussed in section 2.1 the application of our popula- 
tion synthesis model, which assumes a continuos mass function, to a 
real astronomical population is expected to be significantly disturbed by 
population statistics. Whereas for galaxies the mass function is sampled 
sufficiently dense to justify a quasi- analytic treatment this might not be 
the case for OB associations or even worth open clusters. We there- 
fore studied the statistical errors as function of the richness by means 
of Monte Carlo version of the model. Figure 5 shows the spread of the 
resulting interstellar 26 Al masses at time of the emission maximum in 
dependence of the richness. 
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Figure 5 Statistical error 
of the interstellar 26 Al 
mass as function of the 
richness of the association 
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For populations richer than 100 stars initially more massive than 8M Q 
the statistical error drops below the theoretical uncertainties. In ad- 
dition, the Monte Carlo code allows the determination of probability 
density functions for direct application of the association model to ob- 
served populations by means of a Bayesian analysis. 

4. 1.8 MeV from Cygnus 

The Cygnus region is the most significant isolated structure in the 
COMPTEL 1.8 MeV maps beyond the inner galaxy. Figure 6 shows the 
Cygnus OB associations superimposed as circles on the latest COMP- 
TEL Maximum Entropy image of the Cygnus region [22] . 
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As a first test we applied our model to the given OB associations en- 
circled in Figure 6 and generated an 1.809 MeV intensity model, which 
for simplicity neglects the gas dynamic effects discussed in section 2.2. 
The population synthesis is based on a database built from recent liter- 
ature on population studies of the Cygnus region. Figure 7 shows the 
resulting intensity model. At first glance this model already reproduces 
the observed emission peak and the extension to higher longitudes. But, 
significant discrepancies are still remaining. First, the overall flux is un- 
derestimated by the model. After a variation of the slope of the initial 
mass spectrum turned out to be insignificant in resolving this problem, 
the most plausible explanation is an underestimation of populations due 
to strong visual extinction towards Cygnus. Second, the modelled in- 
tensity distribution is significantly narrower than the observed pattern, 
which could be understood if one reminds the neglect of gas dynamic 
effects due to bubble formation. 

5. Summary 

We have presented an extension of population synthesis studies to 
the domain of 7-ray line astronomy. The discussed theoretical uncer- 
tainties are far from being insignificant. In addition, it is questionable 
if the observation and modelling of OB associations as laboratories for 
testing specific aspects of the theoretical understanding of the chemi- 
cal evolution leading to interstellar enrichment of radioactive isotopes 
are sufficient due to possible tremendous statistical uncertainties. For 
stellar populations of massive stars being richer than 100 objects our 
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model shows that theoretical uncertainties begin to dominate. So, one 
can conclude that indeed areas with rich associations may be used for 
testing theoretical aspects. At least, the observations can be used to 
check the consistency by extracting constrains from combined analyses 
at different wavelengths. 

It was shown explicitly that the combined analysis of free-free emission 
and the observed 1.809 MeV intensity may serve as a sensitive age indi- 
cator for young stellar populations. Furthermore, the rather simplistic 
application of the model to the 1.809 MeV emission in the Cygnus region 
already showed the potential in extracting constrains from the 7-ray line 
observations. However, due to the neglect of additional possible source 
candidates such as peculiar massive close binary systems these constrains 
are not very stringent up to now. 
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